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We discuss recent progress of understanding the phenomena observed in high Tc cuprates by study- 
ing the d-wave resonating-valence-bond (d-RVB) based trial wave functions for the 2-dimensional 
t-J-type models. Treat exactly the strong correlation between electrons by numerical approach, we 
examine the evolution of ground states described by different variational wave functions and proper- 
ties of the quasi-hole and -particle excitations of the d-RVB superconducting (SC) state. Properties 
related to the Fermi surface geometry deduced from quasi-hole energy dispersion of the SC state is 
shown to be consistent with the observation by photoemission spectroscopy. With the calculated 
spectral weights (SW's) for adding and removing an electron, we found not only the anti-correlation 
between conductance peak height and width between peaks seen in tunneling experiments, but also 
unique properties due to strong correlation which need to be verified by future experiments. 



The emergence of superconductivity as holes or elec- 
trons doped into the insulating parent compounds is one 
of the fascinating facts seen in high Tc cuprates^. In 
addition to that, there are many other intriguing phe- 
nomena observed in the underdoped regime which may 
indicate new and even exotic phasesi. Phases at dif- 
ferent doping levels have been characterized experimen- 
tally by studying in detail the low-energy, single-particle 
spectra by applying, for example, angle-resolved photoe- 
mission spectroscopy (ARPES) and scanning tunneling 
microscopy/spectroscopy (STM/STS)Si2i^i^. 

To understand these phenomena, we have been study- 
ing the trial wave functions (TWF's), in the form of 
the resonating-valence-bond (RVB) state proposed by 
Andersoni, for the t-J-typc model with different concen- 
tration of doped holes or electrons. Some recent progress 
of this work is discussed in this paper. We start by dis- 
cussing the different variational WF's describing ground 
states at different doping levels. Although our study is 
for only the uniform, 2-dimensional case representing the 
ideal Cu02 plane without any disorders (due to lattice, 
inhomogeneous doping etc.), the phases we found, in par- 
ticular, the antiferromag netic metal (AFMM)2, and the 
coexisting state of superconductivity and antiferromag- 
netism (AFMM-|-SC)^ (see FigUl below)- have recently 
been observed, for the first time, in the Hg based 5-layer 
cupratesiS. 

We have also extended the TWF approach to exam- 
ine properties of the low-energy excitation of the d-wave 
RVB (d-RVB) state in the superconducting regime. In 
the second part of this paper, we present the doping de- 
pendence of Fermi surfaces, extracted from the quasi-hole 
energy dispersion, is consistent with results obtained by 
ARPES. Furthermore, we have calculated the spectral 
weights (SW's) of the quasi-hole and -particle excitations 
and made comparison of the quantities derived from them 
with the observation made by STM/STS. 

Evolution of ground states with doping. Let us 
first discuss the TWF's for ground states of the t- J-type 
Hamiltonian H=J2i j„ —tij c\^Cja + H.c.+Hj with dif- 



ferent numbers of holes or electrons. Here i,j in the 
summation represents nearest-neighbor (n.n.) sites the 
electron hops (with amplitude t), the second n.n. {t'), 
third n.n. site (t") pairs and the Heisenberg exchange 
term -ff/=^X]{i j) (^i ' ^ between n.n. elec- 

trons. Note that Cia- in iJ creates different kinds of holes 
from single-electron-occupied sites: empty holes (Oe-hole) 
for hole doping and two- electron- occupied holes (2e-hole) 
for electron dopingii. Namely, Cia- is actually equal to 
CiCT(l ~ 'T'j,-cr) or Ci^-a-nia foT holc or electron doped case, 
respectively. 

Taking the point of view of doping the Mott insula- 
tor with antiferromagnetic long range order (AF LRO), 
we discuss TWF's determined variationally at different 
doping levels from the undoped case. The TWF's we 
constructed may include three mean- field parameters, de- 
pending on the doping concentration: the staggered mag- 
netization ms={S^)=—{Sg), where the lattice is divided 
into two sublattices when AF LRO is present; the uni- 
form bond order parameters X={J2a ^L'^jo'); and d-RVB 
order A=(cj|Ci| — Cj|Ci|) if i and j are n.n. sites in x- 
direction and — A in y-direction. 

The TWF for lightly doped case is a generalization of 
the single-hole WF with explicit AF LRO first written 
down by Lee and Shih^^. In contrast to states at higher 
doping (see below), WF in this regime is constructed 
solely from the optimized one at half-filling, i.e. with 
Hamiltonian Hj, 

\^o)=P,[ E B^bl/_^,)f^/'\0). (1) 

k[SBZ] 

Here coeflticients A{B)k =[+(-)£^k + Ck]/Ak with 
^k=(Ck,o + A^)^/2, a,o= [4 + (msv)^]^ are energy 
dispersions for the two spin density wave (SDW) bands 
with ek——{cos]ix + cosky). Oko- and bko- represent the 
operators of the lower and upper SDW bands, respec- 
tively, and are related to the original electron opera- 
tors Ck(T and Ck+Qcr with Q=(7r, tt) set for the com- 
mensurate SDW state. Ak=A„(cos fc^; — cosfcj^) for the 
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d-RVB order parameter. The projection operator Pd 
enforces the constraint of one electron per site. At 
half-filhng, Ng- number of lattice sites- equals the to- 
tal number of electrons. Notice that the sum in |^'o) 
is taken over sublattice BZ (SBZ). There are two varia- 
tional parameters: A^, and rusv in this WF. It has been 
shown that this WF produces both the staggered mag- 
netization, {'m)=Ns~^'^j^{—iy^ , and optimized energy 
very similar to those obtained by quantum Monte Carlo 
approach^. 

When a hole is doped or an electron is removed from 
the Mott "vacuum" jvl/o), a pair must be broken with an 
unpaired spin left. Thus it is quite natural to have the 
following Lee-Shih WF for a single doped hole and a lone 
up spin, for example, 

1*1 (q)) = Pd 4 

[k#q] 

where the hole momentum q is excluded from the sum if 
q is within the SBZ, otherwise, q — Q is excluded. Note 
that this WF does not contain any information about 
hoppings of the doped hole or electron. Nevertheless, 
the effect of n.n. hopping is included in the RVB uni- 
form bond which describes the large quantum fluctua- 
tion and spin singlet formation. There is also no need 
to introduce t' and t" in the TWF as they are compat- 
ible with AF LRO. Variational energies calculated vary 
with q and we may obtain the energy dispersion for the 
t-t'-t"-J modelSiiS. The ground state with one hole is 
found to have momentum (7r/2,7r/2) = Q/2. As shown 
in Ref.^, this dispersion relation is still followed when the 
hole number is increased. It is interesting to note that 
states in the hole- and electron-doped cases are in one- 
to-one correspondence after a local particle-hole transfor- 
mation c] _^ is made. Also, because of the Fermi 
statistics, the exchange of a single spin with a 2e-hole has 
an extra minus sign as compared to the Oe-hole. Hence, 
the only difference between the hole and electron doped 
t't'-t"-J models is t' /t -t' /t and t" /t -t" /t after 
we change the Cia on one sublattice sites to —cict-'^. With 
all these, we treat the hole and electron doped cases in 
the same manner. As for the case of having an extra 
up-spin electron with momentum q doped into the half- 
filled state, the energy dispersion can be calculated with 
this same l^*!). The ground state is here at momen- 
tum (tt, 0). This result agrees well with that obtained by 
self-consistent Born approximation (see the discussion in 
Ref.^). Dispersions for single- hole and -electron cases 
turn out to be simply the combination of the mean-field 
band at half-filling and the coherent hopping o^i^^ . Using 
l^*!), we calculated the momentum distribution function 
(MDF) nk — (cj^^'^kcr) for the ground state of, say, a sin- 
gle Oe-hole with momentum Q/2. We observed the same 
dips or pockets at Q/2 and anti-dips at -Q/2 found for 
the exact results of 32 sites (see Ref.^). 

We may generalize the WF to describe states with two 
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FIG. 1: Ground state phase diagram for the t-t'-t"-J model 
with {t' ,t")/t—{— 0.2,0.1). Shaded and crossed-line regions 
represent schematic doping regions of AFMM and coexisting 
AFMM-I-SC states, respectively, (m) represents the staggered 
magnetization value (vertical axis on the left) and P^""' ( axis 
on the right) for SC pairing amplitude (see text for their defi- 
nition). The results are obtained on the lattice of size 12 x 12. 



doped holes by taking out the unpaired spin from l^*!). 
The lowest energy state turns out to be the one with zero 
total momentum and Sz=0. Namely, the TWF for two 
holes with momenta q and — q is 

|*2)=Pd[E'(^kat^«^-ki + 5k&l,T&^ki)]'"'^/'^"'|0)- 

k 

.(2) 

Note that momentum q is not included in the summation. 
It is most surprising to find that, although \'i'2) has zero 
total momentum irrespective of q, its energy varies with 
the missing momentum or the hole momentum q. Same 
type of TWF for four holes with momenta ±q and ±q' 
both excluded from the summation can be constructed. 
The energy dispersion for two holes doped into the half- 
filled state is also almost identical to that of a single 
hole and the minimum is at Q/Qm. The MDF for this 
state has both dips at Q/2 and -Q/2. This is in good 
agreement with the exact result for the t-t'-t"-J model 
with 2 holes in 32 sites. For 2 electron doped case, the 
dispersion turns out to be very similar to that of a single 
doped electron^. The state with momentum (tt, 0) has 
the lowest energy for two electrons. 

Our results reproduce the contrasting behav- 
iors between lightly hole and electron doped [e.g. 
Nd2-xCe.j;Cu02 (NCCO)] high cuprates revealed by 
high- resolution ARPES^'^'"'. While a small hole patch is 
observed to be at the center in the nodal direction of 
the BZ [i.e. (Q/2)] in hole doped systems, small electron 
patches centered at (7r,0) and (0,7r) are observed in 
lightly doped NCCO. 

Focusing on the hole doped case, we found that 
this type of TWF with even number of holes have 
AF LRO but very little superconducting pairing 
correlations- namely, P^^'^, the averaged value of the 
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long-range part of d-wave pair-pair correlation func- 
tion, Pd(r) = (iV,)-l J2^{^l^^+r) wllCrC Ai=C,,| (c,+x,i + 



for Irl > 2 in the 12 x 12 lattice is 



only about 10"'* (x and y are unite vectors in the x and y 
directions, respectively). The tiny pairing strength is also 
indicated by the very small value of n.n. hole-hole cor- 
relation function, (Ns)~^ ~ ~ 'T-i+r)) with rii 
the electron number at site i, as shown in Fig. Efd). This 
state with large staggered magnetization, small Fermi 
patch, but no SC pairing behaves like an AFMM [see 
Fig. Ill and Fig. |3a) and (d)]^. Note that this state has 
been found in the inner planes of the Hg-based 5-layered 
cuprates H gBa2Ca4^Cu^Oi2+r^- 

As more holes are doped into the system, the ground 
state may be switched from AFMM to the one described 
by WF, 
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(3) 
and 



of electrons 
y) containing 3 new 



where iVg is now the number 

±^k = ±Ck,l= ±[(c0sk:r -|-C0skj^)2 

COS ka; cos ky — (cos -f cos 2k 
variational parameters /i^, for chemical potential, and 
for long range hoppings. This state, which we 
denote as AFMM-I-SC, has stronger superconductivity 
than that of AFMM state [P^°""= is about 5 x IQ-^, see 
Fig. ^ and finite correlation between n.n. holes, see Fig. 
Eld)] and more extended MDF [Fig. Elb)]. SC and AF 
ordered states are uniformly mixed or coexisted in the 
AFMM-fSC. 

Further increasing the number of doped holes in the 
system, the ground state becomes the RVB-SC state, 
namely, 
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FIG. 2: Plot of MDF's for 12 x 12 sites obtained by (a) 
AFMM (j;=0.014), (b) AFMM-hSC (a;=0.042), and (c) RVB- 
SC (x=0.069) states with the same parameters as that in ob- 
taining results shown in Fig.l. (d) The hole-hole correlation 
functions for different WF's, as denoted by different symbols. 
Vertical axis on the right (for AFMM) is with a finer scale 
than that on the left. 



the effects of strong correlation exactly by investigat- 
ing the SW's of the uniform SC state on finite square 
lattices*^. 

We consider the quasi-particle excitation defined by 
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where Ok = fk/^k — [E^ — ek)/Ak with Vk and Uk d- 
wave SC coherent factors. Here, ek — ^k,i{'msv — 0) in 
the unit of n.n. hopping amplitude. As shown in Fig. ^ 
and Fig. 12c) and (d), this state shows even stronger SC 
pairing amplitude and has a large Fermi surface^. 

The evolution of states, from AF LRO, AFMM, 
AFMM-fSC to RVB-SC, we described above for ideal 
2-dimensional system have been realized recently in mul- 
tilayer system which is considered to have perfectly flat 
Cu02 plane&lSi. Indeed, the ground state phase dia- 
gram. Fig. ^ as determined for proper set of parameters 
(t', t")lt is similar to the one by experiments. 

Properties of the quasi-hole and -particle ex- 
citations. We then discuss our attempt to apply the 
TWF approach to study the properties of quasi-hole and 
-particle excitations of the RVB-SC ground state. We 
have calculated the quasi-hole excitation spectrum of the 
SC state and deduced the normal-state Fermi surface. 
Without any disorder in the system, we also examined 
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By applying the t-t'-t"-J Hamiltonian to \Nf. — 1), we 
calculate its excitation energy for each momentum and 
then fit the results by the excitation dispersion, i?k= 
{el + A2)i/2 .^vith = Ck,i(wa„ = 0). The Fermi 
surface is then determined by the zero energy contour 
of ek. The doping dependence of the hole density, de- 
duced from the area of the determined Fermi surface 
is shown in Fig. Of a). With different sets of {t',t")/t 
which are believed to be material dependent- for exam- 
ple, in Fig. Ela), (-0.3,0.2) for Ca2-xNa^Cu02Cl2 and 
(-0.1,0.05) for La2-xSrxCu02i^Li^, the results are con- 
sistent with that determined by ARPES. We may also 
examine the Fermi surface position along the nodal direc- 
tion, i.e. along (0,0)-(7r, tt), at each hole doping density 
[see Fig. E^b)]. As can be seen clearly, the Fermi sur- 
face position moves more prominently for the case with 
larger {t',t")/t values, again consistent with the experi- 
mental observation (see also Refi^ for Bi-based bilayer 
cuprates) . 

Let us then move to calculate the SW for adding (and 
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FIG. 3: Doping dependence of (a) the hole density xfs and 
(b) the position at the nodal point obtained from the 
fitted quasi-hole dispersion. (See text for the detail.) Solid 
hne in (a) represents the relation xfs ~ x. Solid circles and 
triangles are numerical results for {t' ,t") /t={-0.3,0.2) and (- 
0.1,0.05), respectively. Dashed and solid lines in (b) are linear 
fitting to these two sets of results. Empty symbols represent 
experimental data extracted from ARPES, with Na — CCOC 
for Ca2-^Na^Cu02Cl2 and LSCO for La2-^Sr^Cu02i^. 

removing) one electron defined by 

^+(-) ^ l(^e + (-)l|cL(c,..)|iVe) ? 

Applying identities for projection operator, 

[Cka,Pd]Pd = 0; 

PdCwAcl,,Pd] - Pd[^5]e'^'''-''^-^'"n,,_,]Pd(8) 

I 

with Rio- the position vector of the i-th spin a in the 
lattice of size Ns and = c\^Cia, we can relate Z^^ 
exactly to the hole doping and MDF as 

ZL = ^^-r^^.. (9) 

On the other hand, reminiscent of what have been 
argued previously by analytic approach^i, we recognize 
that the strong correlation effects becomes apparent only 
in Z^^ at low doping. The effects due to strong correla- 
tion are examined by comparing the coherent SW aver- 
aged over all momenta, i.e. Z~^^ = Z^^/Ng, and the 
incoherent part defined by the relation 



obtained by exact treatment of the projection and by 
using renormalized mean-field theory (RMFT)?i. Here 
riave = X]k"kcr/^s IS the average MDF which should 
always be equal to the electron density of the system. 

The exact results for the 12 x 12 lattice and that 
by RMFT are shown in Fig0Ja). The coherent part 
of Z^^ by RMFT is gtv^ with renormalization factor 
gt — 2a;/(l -|- a;). Completing the momentum sum for the 
coherence factor, the average result is x{\ — x)/{\ + x) 
and, thus, n^^^l = [(1 - a;)2]/2(l + x)^. The RMFT re- 
sults are plotted in Fig^a) (dashed and dotted Hues) in 
comparison with the exact ones. As is shown there, while 
the numerical riave is indeed equal to the electron den- 
sity, the exact incoherent SW for removing an electron 
is less than the RMFT result. The difference becomes 
more significant as hole doping level is reduced. Inter- 
estingly, this behavior is independent of the (t', t")/t val- 
ues [represented by solid and empty symbols in Fig^Ja)] 
which correspond to very different doping dependence of 
the Fermi surface shape and also the density of states. 
By contrast, the average values of Z^^ calculated exactly 
(not shown) and by RMFT are identical due to Eq.Q. 

As a brief digression to electron doped case, it is 
straightforward to show that, applying the hole-particle 
transformation to Eq. , SW for removing an electron in 
electron doped system satisfies the relation .^kCT='^kcr — 
(1 — x)/2 rigorously. The doping dependences of Z^^^ 
and Uave are shown in Fig. 4(b). The numerical results 
is exactly equal to RMFT ones indicating that there is 
no incoherent part of SW for removing an electron in the 
electron doped system. 

To make a comparison with tunneling experiments, we 
then concentrate on the SW's as a function of the excited- 
state energy. To reduce the effect of finite size, we define 
the sum of Z^^/N , over momentum k which has energy 
within E - AE/2 and E + AE/2, as g{E) which could 
be viewed, approximately, as proportional to the conduc- 
tance at low energy E. We plot g{E) in Fig. up to 
about the energy where peaks appear for lattices of size 
12 X 12 with AE/t =0.3 for various dopings. (The energy 
interval we chose is to reduce the effects due to lattice size 
and also to include about the same number of momen- 
tum points below and above zero energy.) To make sure 
our treatment is correct, we have also applied the same 
analysis to the d-wave BCS (d-BCS) state and found the 
ideal BCS result is hardly distorted by the finite size. 
Therefore, with the reasonable finite-size dependence, we 
obtain indeed the F-shape d-wave gap near zero energy. 
The width between peak positions is also roughly equal 
to two times of the gap value deduced from the excita- 
tion energy. Looking at the results closely, while g{E) 
may indeed be about the same at the opposite sides in 
the very vicinity of zero energy as suggested in ReffSi, 
g{E) for removing an electron is always larger than that 
for adding an electron at higher energy near that of the 
peak. With decreased doping, the ratio of g{E) at neg- 
ative and positive energies enhances quite dramatically, 
e.g. from a;=0.125 to 0.056, g{—A)/g{A) at the corre- 
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FIG. 4: The doping dependence of average SW and MDF 
of RVB-SC state in hole doped (a) and electron doped (b) 
systems, obtained numerically for the 12 x 12 lattice and by 
RMFT. Numerical and expected Uave. are denoted by the cir- 
cles and the solid line, respectively. The squares (triangles), 
connected by dashed (dotted) line as the guide for the eye, are 
for exact results of Z^„e \n^ave only for hole doped case, ex- 
tracted using Eg. nun with riaue calculated numerically]. The 
dashed and dotted lines without data points represent results 
by RMFT. All solid symbols in (a) are resuhs for {t' ,t")/t={- 
0.3,0.2) and empty ones for (0,0). Results in (b) are all for 
(t',t")A=(0.3,-0.2). 



sponding peak energy A (in units of t) increases from 
1.96 to 2.73. Similar behaviors are found for the case 
with vanishing {t',t")/t, as shown in the inset of Fig. |S1 
In contrast to this, for the d-BCS case in the same finite 
lattices there is almost no change of the ratio within the 
gapi&. The numerical results thus tells us features due 
to strong correlation, i.e. the particle- hole asymmetry 
of average conductance exists even within the gap region 
and gets enhanced with underdoping, which are not yet 
fully explored in the tunneling experiment. 

We can also see, from Fig. |S1 the correlation between 
heights of the spectral weight peak and the gap size as 
doping level is varied. For the two doping level shown 
there, the peak height scales with the pairing amplitude 
but anti-correlates with the gap size. This is in clear 
contrast to the BCS case in which the peak height, pro- 
portional to the SC coherence, scales with the gap size 
as more holes doped into the system. Our result agrees 
qualitatively with what has been extracted from STS 
experiments'*. 

To summarize, we discussed d-KVB based TWF's de- 
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FIG. 5: g{E) for RVB-SC state versus excitation energy E for 
2 different dopings. Results shown here are for {t' ,t")/t={- 
0.3,0.2) and lattice size 12 x 12. The associated excitation gap 
positions are marked by arrows (see text). Inset: Same plots 
for {t' ,t")/t— (0,0). Results for two doping concentration are 
denoted by the same symbols as that in the main figure. 



states have been observed recently in multilayer cuprates. 
We have also determined the doping dependence of both 
the Fermi surface area and Fermi momentum along the 
nodal direction using the fitted quasi-hole dispersion. 
With proper parameters for different real materials, the 
results are consistent with the same quantities observed 
in various high Tc cuprates. By studying SW's of the 
quasi-hole and -particle excitations of the SC state, we 
found, even without any disorder, features qualitatively 
consistent with what has been seen in tunneling experi- 
ments. An exact relation is obtained for the doping de- 
pendence of the spectral weight of removing an electron 
from an electron-doped cuprates. This could be tested 
by ARPES. There are also other results which have to 
be confirmed by future experiments. 
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